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A PRELIMINARY COLLATION OF THE 
THERMODYNAMIC AND TRANSPORT PROPERTIES OF CESIUM 


H. W. Hoffman and T. T. Robin, Jr. 


ABSTRACT 

■ ^ 

This memorandum constitutes a preliminary collation of the thermo¬ 
dynamic and transport properties of saturated liquid and vapor cesium. 

It is intended to provide a convenient and consistent source of prop¬ 
erties for those involved in the design and evaluation of power systems 
based on thermodynamic cycles involving cesium as the working fluid. 

Included in this compilation are the thermodynamic properties — 
physical constants and critical state conditions, sonic velocities, PVT 
relationships (vapor pressure, density, and compressibility), enthalpy, 
heat capacity,- latent heat of vaporization, and entropy, and surface 
tension — and the transport properties — viscosity, thermal conductivity, 
and electrical resistivity. A preliminary evaluation of the data avail¬ 
able indicates the need for additional measurements of both the thermo¬ 
dynamic and transport properties. Thus, for liquid enthalpy, the dif¬ 
ferences between the tabulated relations, while small, are sufficient to 
affect significantly thermodynamic cycle calculations. To a lesser extent, 
further data on the PVT characteristics would be useful, since these data 
are the base for the calculation of vapor enthalpies; only Ewing et al. 
predict vapor pressures for temperatures above 1950 °F. Achener alone- 
reports measurements of the latent heat of vaporization; the remaining 
’'investigators M derive this quantity, with only Ewing showing reasonable 
agreement with the experimental values. The surface tension has been 
obtained only by Tepper to date; experiments by others are in progress, 
and evaluation of need for further studies must await the outcome of these 
current investigations. In respect to the transport properties, serious 
disagreements exist for vapor viscosities and thermal conductivities; and 
the need for additional measurements is evident. Some variation, though 
not as significant, exists for the liquid viscosities and conductivities. 

The results of Ewing et al. are currently recommended for the thermo¬ 
dynamic . properties in that these constitute an internally consistent set. 

A later revision of this memorandum will include a more detailed evalua¬ 
tion of this data set in relation to similar consistent sets by Heimel 
and by Achener et al. Values for superheated vapor and the nuclear 
properties of cesium will also be included in the final collection. 
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IRTRODUCTIOR 

The alkali liquid metals have been proposed as the coolant and 
thermodynamic medium in space power applications utilizing nuclear 
reactor or isotopic heat sources. These fluids are also of interest in 
the fast-breeder reactor program where they serve as the core coolant. 

The work summarized herein is a part of Phase I of an analytical com¬ 
parison being carried out at the Oak Ridge Rational Laboratory for the 
Rational Aeronautics and Space Administration under AEG Interagency 
Agreement 40 - 98 - 66 , RASA Order W- 12 , 353 > on "the relative merits of cesium 
and potassium as working fluids for Rankine-cycle space power plants. 

The thermal properties of the alkali metals are necessary both for 
cycle analyses and for the correlation and interpretation of experimen¬ 
tal heat-transfer data. This report compiles the thermodynamic and 

% 

transport properties of only one of these fluids - cesium; a subsequent 
addition will cover potassium. A preliminary evaluation of the available 
data is included. A final report issued at the conclusion of this pro¬ 
gram will contain (l) a set of recommended equations,‘ and ( 2 ) an indica¬ 
tion of areas in which further experimental measurements are necessary. 

The thermodynamic and transport properties of saturated liquid and 
vapor cesium are presented in the sections which follow. The data are 
given in two forms: (l) graphs of the specific property over the temper¬ 
ature span 1100 to 2500 °R constituting a preliminary set of values recom¬ 
mended for design calculations; and ( 2 ) tabulations of available equations 
for each property including temperature range, reported error, indication 
of the method of measurement, the investigator and the date of the investi¬ 
gation, and the reference. Each of the two data sections is preceded by a 
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brief discussion of the properties within- that section and a preliminary 
evaluation and justification of the recommended data. These discussions 

include — particularly in the section on transport properties — comments 

• ♦ ■■ 

on properties which are not given in ensuing data tabulation,; the primary 

i 

purpose for this has been to provide a reasonably complete perspective on 
cesium properties. Further, little detail is included in these comments; 
for further elaboration-reference should be made to the original sources. 
While the application of this collection is in the design of nuclear 

systems utilizing cesium as the reactor coolant.and/or thermodynamic 

* \ _ 

medium, no information on; nuclear properties is included,. The possibility 

of such data appearing in a final version of this report will be considered. 

For the thermodynamic properties, the results of Ewing et al . 1 have 
been recommended; the basis for this choice is discussed further in the 
section immediately following. Subsequent to the initial writing of this 
memorandum, Achener et al . 35 issued a report in which, using their meas¬ 
urements as a basis, a self-consistent set of thermodynamic properties was 
also developed. A brief comparison of the two property sets indicates 
some differences. The significance of these differences has not yet been 
evaluated but will be included in the final version of this memorandum* 

• % i 

Superheat values for cesium vapor are not included in this initial 
cesium properties summary. These should be obtained by reference to 
Appendix B of Ref. 1 ; for the convenience of users of this data collection, 
tabulations or plots of properties of superheated cesium vapor may be 
appended to the final report. 
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THERMODYMMIC PROPERTIES 

Introduction 

The thermodynamic properties characterize a substance with respect 
to equilibrium states. In this category, we include the critical state 
conditions and sonic velocity, the PVT relationships (vapor pressure, 
density, compressibility), the enthalpy, heat capacity, latent heats of 

i 

phase change, and entropy, and the surface tension. 

Ewing and co-workers 1 and Heimel 15 have both generated internally 
consistent sets for the thermodynamic properties employing approximately 
the same procedures. The former list saturation values between 1250 and 
2350 °F (limiting their analysis to pressures between about 1 atm and 
3 ^ atm), while the latter tabulates saturation properties over the span 
— 280 °F to 2240 °F. In the common range, there are some generally small 
differences between the two data sets. Ewing et al. have also calculated 
superheat values (Appendix B, Ref . l) for cesium vapor between 1250 and 
2550 °F (to a maximum pressure of 3^ atm) using as the starting point for 
a particular property the absolute value of that property for either the 
saturated liquid or the monomeric gas. In general, the agreement between 
values calculated by either path was quite good (see Ref. 1 , p. 21 ). 

This is discussed further in a following section on enthalpy and entropy. 

Because of the wider range, the Ewing et al. results are currently recom- 

* * 

mended. 

One further comment is necessary. It is known that the vapors of 
many metals do not behave as normal gases because the molecular weights of 
such vapors change as conditions change. This is generally believed to 
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arise from variations in the relative amounts of monomer, dimer, trimer, 
etc. species in the vapor. Since a certain amount of energy is tied up 
in these polymeric states, it is necessary to account for these vapor 
species composition effects in calculating the thermodynamic properties 
of the vapor. ' When this is done, the properties are known as equilibrium 
properties . On the other hand, when processes occur rapidly, equilibrium 
between the several vapor species may not be achieved; e.g., as in the 
nozzle expansion of a vapor. Calculations are then made on the basis of a 
fixed vapor composition and the.properties are known as frozen properties . 
Generally, under conditions of shifting equilibrium, bracketing calculations 

are made using the extremes of equilibrium properties and frozen properties. 

* * * 

Physical Constants and Critical Properties 

A brief tabulation of some of the physical constants for cesium are 

given in Fig. 1; nuclear properties are omitted but will‘.be included later. 

The critical properties of a. fluid (T , P , and v ) are important in the 

c c c 

design of vapor compressors and turbines and useful in the estimation of 
the transport properties. Further, some generalization of the properties 

V 

of different substances is possible through the use of the reduced state 
concept (ratio of the existing state to the critical state). 

Wo experimental values are available for the critical properties of 
cesium; all values listed in Table 1 were estimated through a variety of 

i 

methods. As a result, no recommendation as to best values is given. 

Critical temperatures range from 1838 to 2157 °K (the extremes of the esti- 

* 1 

mates by Tepper et al. 17 ) with a mean of 2055 °K, critical pressures vary 
from 140 to 170 atm, and critical volumes from 0.311 to O.38O liter/g-mole. 
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Sonic Velocities 

The equilibrium and frozen sonic velocities for saturated vapor are 
shown in Fig. 2 . The curves plotted were taken from the calculations of 
Shapiro and Meisl 31 as reported by Weatherford et al . 26 Values based on 
experimental measurements are not available. 

PVT Relationships 

Vapor Pressure . Six vapor pressure equations are shown in Table 2 . 
Some,early data has been omitted due to unavailability of the • references 
at this time; these may be included in the final version of this memorandum. 
The recommended equation, that of Ewing et al ., 1 is shown as Eq. ( 3 ) in - 
Table 2 and plotted in Fig. 3 * 

The six tabulated vapor pressure equations show the following com¬ 
parison at the three listed temperatures: 


Vapor Pressure 


Equation 


(psia) 


850 °F 

1700 °F ' 

2300°F 

(i) 

1.01 

99-5 

462.7 

( 2 ) 

0.98 

100.4 

472.2 

( 3 ) 

1 - 03 

3 . 00.8 

452.2 

( 4 ) 

1.05 

105-3 

^ 93-7 

( 5 ) 

0.84 

65-3 

252.9 

( 6 ). 

1.05 

103.6 

449.1 


The relation of Taylor and Langmuir 14 [Eq. ( 5 )] deviates significantly 
over the entire temperature range from the other five expressions; 
Bonillas 13 results are somewhat higher. Tepper , 12 Achener , 2 and Ewing 1 
are in close agreement except at 2300 °F which exceeds the upper temperature 
limit of the data of the first two. While the lower limit of Ewing 1 s 1 data 
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is 1236 °F as compared to 850 °F for Tepper, 12 893 °E for Achener, 754 °F for 
Bonilla, and 85°F for Heimel, 15 the agreement in the calculated vapor 
pressure remains reasonable. 

To determine the effect of vapor pressure on thermodynamic cycle 
calculations, the equations of Tepper and of Ewing were used in an ideal 
RanHine-cycle calculation. Turbine inlet conditions were assumed as 
2250 °F and 293*92 psia; the exhaust condition was 58-78 psia. Using 
Eq. (l), the cycle efficiency was found to be 20 . 007 1 ° with an exhaust 
quality of 0 . 86897 - Correspondingly, Eq. (3) gave a cycle efficiency of 
19.981$.with O.8708O exhaust quality. Thus, the effect of the particular 
vapor pressure equation elected (using relationships in fair agreement) 
on the engineering cycle calculations does not appear large. 

Density . For liquid cesium, the equation of Ewing et al. 1 [Table 3 > 
Eq. (3); and Fig. 4 ] is recommended. This equation represents the best 
fit to the experimental data obtained by Ewing, Achener, 2 ' 3 and Tepper 
et al. 5 over the temperature range from the melting point to 2300 °F. 

For cesium vapor, Ewing et al.f [Table 4 , Eq. (l); and Fig. 5] is 
recommended. The virial equation of Ewing predicts experimental data to 
within ± 0 . 43 $; the empirical equation of Tepper et al . 17 correlates their 
data with a maximum deviation of 7 - 99 $- Saturation values only are plotted 
in Fig. 5 • 

Compressibility . The compressibility, z = pV/RT, is often a useful 
quantity in calculations. The calculated values of Ewing et al. 1 arc 
again recommended, z is plotted in Fig. 6 for cesium vapor at saturation 
conditions only; values for superheated vapor are to be found in Appendix B 


of Ref. 1 . 
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Enthalpy and Entropy- 

Enthalpy . The five equations shown for the liquid enthalpy in Table 5 
agree reasonably with respect to enthalpy differences over the range 1000 
to 2500 °R. The Ewing et al. 1 relationship [Table 5 , Eq. ( 4 ); and Fig. 7 ] 
is recommended at this time. 

Equation ( 4 ) is based on the absolute properties of solid cesium at 
77 . 0 °F as measured by Hultgren et al. 28 and the data of Lemmon et al. 4 
for the latent heat of fusion. An anomoly in the heat-content curve for 
the liquid exists between 77 and 570 °F, and the enthalpy change as obtained 
by Lemmon et al. was used for this region. A constant value for the 
specific heat of the liquid of 0 . 05^83 Btu/lb m * °F was used above‘ 570 °F 
(see discussion on c ), The values tabulated by Ewing (Appendixes A and 
B, Ref. l) were based on monomeric gas path calculations — as discussed 
below — rather than on Eq. ( 4 ) by subtracting at each temperature the 
enthalpy of vaporization. The resulting differences in enthalpy were 
small (0.1 to 2.0$). 

The enthalpy of saturated vapor is plotted in Fig. 8. Ewing et al. 1 

obtained the vapor enthalpy by computation along isotherms using their > 

« 

« 

PVT equation. Thus, for the monomeric gas. 
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where is the molecular weight of monomer (132.91) and the remaining 
• « 

qunatities are as shown in Table 4 , Eq. (l). The reference value for the 
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monomeric gas at 1 atm (relative to the - solid crystal at 0°R) is given by: 

H° = 252.18 + 0.037361 T + 2k0 e -31 } 290/l 


where H is in units of Btu/lb and T. in °R. This equation is based on 


the standard properties of Evans et al. 33 over the temperature range 

0 to 3100°F and on the enthalpy of vaporization to 0°R (18.62 mean kcal/ 

mole) as derived by Ewing. 1 It is noted that the value for used by 

S 


Ewing differs somewhat from that employed by Heimel. 
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Latent Heat of Vaporization . The equation of Ewing et al. 1 is 
recommended; this is shown as Eq. (2) in. Table 6 and is plotted in Fig. 9* 
This relation was derived by differentiation of the vapor pressure equa¬ 
tion [Eq. (3), Table 2] and substitution into the Clausius-Clapeyron 
equation. A value of v at each temperature was obtained from Eq. (3), 

Jo 

Table 3> and of v from the virial equation of state [Eq. (l), Table 4]. 

s 

Of the four equations presented, only that of Achener 2 is based on 
experimental measurements. The other three equations — Ewing et al., 1 
Tepper et al.f 1T and Dillon 16 — obtain the latent heat of vaporization 
by the procedure discussed above. Only Ewing shows good agreement (l to 


3$) with the Achener data. 

. Entropy . • Again, Ewing et al. 1 is recommended [Table 7* Eq. (2); and 
Fig. 10]. As discussed with respect to the enthalpythe values tabulated 
by Ewing for the entropy of the saturated liquid were based on the prop¬ 
erties of the monomeric gas at 1 atm; these are the values plotted in 
Fig. 9. Equation (2) was used in comparison calculations with the prop¬ 
erties of the saturated liquid as the starting point; note that this 
entropy is given relative to the solid at 0°R. The resulting differences 
in-entropy were email (0 to 0.8$). 
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The discussion on the vapor entropy follows that given above for the 

t 

vapor enthalpy. Thus, for the momomeric gas. 



whe re 

S° = 0.037361 &zT + 0.080604 + 0.371 e ’ 2 9;>598/T , 

§ 

for p in atm, T in °R, V in ft 3 /lh-mole, and S in Btu/lb m *°F. Vapor 
enthalpies tabulated by Ewing et al. (Appendix A, Ref. l) for saturated 
vapor are plotted in Fig. 11 . 

The entropy of vaporization at each saturation point is obtained by 
dividing the appropriate enthalpy of vaporization by the absolute temper¬ 
ature . 

Specific Heat . Six equations are given in Table 8 for the specific 
heat of saturated liquid. With the exception of Ewing et al. 1 [Eq. ( 5 )]/ 
these were obtained by differentiation of derived heat content equations. 
Ewing f s value of O.O5683 Btu/lbxn*°F was obtained through a trial-and-error 
procedure to provide internal third law consistency with their PVT data. 
Note that this value lies intermediate to those of Tepper et al. 12 and 
Lemmon et al. 4 and. further, that an anomaly exists in c in the tem- 
perature range between 212 and 572 6 F (see Lemmon et al.. Ref. 4 , p. 12 ). 
Ewing et al. is recommended (Fig. 12 ). 

The specific heat of cesium vapor (Table 9) was determined by Ewing 
et al . 1 following the procedure outlined above for enthalpy and entropy. 
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For this property, however, values based on the saturated liquid were 
5*3 to 31$ lower than those calculated by the monomeric gas path. The 
latter values are plotted in Fig. 13* Shapiro and Meisl 31 report the 
frozen specific heat of saturated vapor to be 0.0371 Btu/lb *°F* 

To evaluate the effect of liquid c; on cycle efficiencies, an ideal 

xv 

Rankine cycle was assumed; and calculations were made for two different 
values of c „ — O.O5683 Btu lb • °F (as recommended by Ewing et al. 1 ) and 

p,j0 m x 

0.0600 Btu/lb ••°F (as reported by Lemmon et al. 4 ). The Ewing relations 
were used for all vapor conditions. The turbine inlet conditions were 
taken to be 2250°F and 293*92 psia; the turbine exhaust was at 58*78 psia 


For the first case, the cycle efficiency was determined to be 19*981$ with 
a turbine exhaust quality of 0.870800; the second case yielded an effi¬ 
ciency of 24.011$ and.an exhaust quality of O.829166. These differences 
are sufficiently significant to suggest the need for further experimental 


measurements on the enthalpy of liquid cesium. 


Surface Tension 

The interfacial tensions (solid-liquid and vapor-liquid) of the alkali 
liquid metals are of importance in situations where free liquid surfaces 
exist, where liquid transfer is by capillary action (heat pipes), and 
where heat transfer occurs by boiling. The only currently available data 
for cesium are those of Tepper et al. 17 as shown in Table 10 and in 
Fig. 14. Achener 34 is currently studying lithium and sodium and plans 
measurements with cesium. Cooke (OREL) has also started experiments with 
cesium to obtain data to about 2000°F. All three investigators are using 
the maximum bubble pressure technique.* : 





Fig. 1. Cesium: Physical Constants 

* 


Chemical Formula: Cs 

Molecular Weight: - 132. 91 

Melting Temperature: 83.-l°F 
Boiling Temperature: 1236°F (l atm) 
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Table 1. Critical Properties 


Equation 

Eq. 

No. 

Temperature 
Range (°F) 

Standard 

Deviation 

Method 

Comments 

Investigator 

(Date) 

Refs? 

Temper ature (°K) 








T c = 2057 



±4o°k 


Error based on cata at l650°C 

Dillon et al. 

(1965) 

16,27 

T = 2000 
c 






Morris (1964) 

7 

T = 2150 
c 






Grosse (i960) 

8 

T = 1996 
c 



• 



Hochman, Bonilla 

(1964) 

18 

T c = 1838 to 2157 





Five methods used for estimating 

Tepper et al. 

(1965) 

17 

Pressi re (atm) 








P = 145 
c 





Equation of Sawhney and Bonilla 

using T of Dillon et al. 
c 

Dillon et al. 

(1965) 

16,27 





* 


Sawhney, Bonilla 

(1962) 

13 

P = l4o 

c 






Morris (1964) 

7 

P = 170 
c 






Grosse (i960) 

8 

Volume (liter/g-mole) 








v = 0.311 

c 



i0.008 



Dillon et al. 

(1965) 

16,27 

v = 0.380 

c 






Morris (1964) 

7 

v c = 0.316 






Grosse (i960) 

8 


* 

See Table of References. 
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Fig. 2. Cesium: Vapor Sonic Velocity. 


Equilibrium Sonic Velocity (ft/sec 





























































































































































































































































































































































































































































































































































































































































































Vapor Pressure (psia) 
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Fig. 3. Cesium: Vapor Pressure. 











































































































































































































































































































































































































































































Table 2. Vaoor Pressure 


Equation 

Ec. 

Nc. 

Temperature 
Range (°F) 

Standard 

Deviation 

Method 

Comments 

Investigator 

(Date) 

Refs* 

logic p = - 663 y — + 3-9011 
(P = atm abs; T = °R) 

( 1 ) 

850-1950 

0 . 72 $ 

Boiling 

No mention of super-heat 
phenomenon or "bumping"; 
for Rb added Zr-chips to 
prevent bumping 

F. Tepper et al. 

( 1963 ) 

12 

logic P = - 668 °j l85 + 5.09475 
(P = psi abs; T = °R) 

( 2 ) 

893-1600 

0 . 68 $ 

Boiling 

No mention of super-heat 
phenomenon 

P. Y. Achener 
(1964) 

2 

logic P = - 7 + 5.87303 - 0.53290 logic T 

(P = atm; T = °R) 

( 3 ) 

1236-2346 

± 0 . 35 # 

Null Point 


C. T. Ewing et al. 

( 1965 ) 

1 

logic P = - 666 y ° 2 + 5.10876 

(P = psi abs; T = °R) 

(4) 

75^1788 

1.44$ 

Boiling 


C. F. Bonilla et al. 

( 1962 ) 

13 

hoiil 

logic p = - -“r + 11.0531 - 1.35 logic T 

(P = mm Kg; T = °K) 

(5) 

84-1340 

300-350 °K: lfo 
350-600 °K: y/o 
6CO—1000°K: % 

Positive Ion 


J. P. Taylor and 

I. Langmuir 

( 1937 ) 

l4 

logic P = - 30 + 7.04453 - 0.915282 logio T 

(P = atm: T = °K) 

(6) 

85-1938 

0 . 5 $ 

Thermodynamic 

Calculation 

Based on zero-degree heat 
of condensation of the 

monomer• 

S. Heimel 

( 1965 ) 

15 


* 

See Table of References. 
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Fig. 4. Cesium: Liquid Density. 








































































































































































































































































































































































































































































































































































































Table 3» Density (Liquid) 


Equation 

Eq. 

No. 

Temperature 
Range (°F) 

Standard 

Deviation 

Method 

Comments 

Investigator 

(Date) 

* 

Refs. 

P = i. 868926 - 0.29211 ( ^) 

- 0-0154436 ( 1q £ 0 ) 

[p = g/cc; T = °F] 

(1) 

100-1950 

0.00236 

g/cc 

Dilatometer 
(constant wt) 


Tepper et al. 
( 1964 ) 

5 

p = 1 . 84 /[l + 1.1755 X 10 4 (T - 82 . 4 ) 

+ 7.656 x 10 - 8 (T - 82. 4 ) 2 ] 

[p = g/cc; T = °F] 

( 2 ) 

83-1671 

0.67$ 

(0.010 gm/cc) 

Pycnometer 
(constant vol) 


Achener ( 1964 ) 

2,3 

p = 124 .181 - 1 .5970 X 1 C ' 2 T 

- 1.6855 X 10 - 6 T 2 

[p = lb m /ft 3 ; T = °R) 

(3) 

83—23OO 

* 0 . 6 % 

Pycnometer 

Equation is best fit to 
combined NRL, AGN, and 

MSAR data. 

Ewing et al. 

(1965) 

1 

Data tabulated as function of temperature 

* 


1112-3218 


Irradiation 

Count of Y activity from 
known isotopes irradi¬ 
ated by thermal neutrons 

Dillon et al. 

(1965) 

16,27 


See Table of References. 
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Fig. 5. Cesium: Density, Saturated Vapor 
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Table 4. Vapor Density 


Equation 


Eq. 

No. 


Temperature 
Range (°F; 


Standard 

Deviation 


Method 


Comments 


Investigator 

(Date) 


Refs. 


PV _ B , C D E 
RT V V 2 V 3 V 4 


logic B = - 3.6200 + 


B < 0 


4000 


+ logi o T 


log C = 3.3551 - ”31^ + 10,^, .x 10 6 


c > 0 


log Id I = 4.1856 + 


83o 


D < 0 


E = +600,000 

[P = atm abs; T = °R; V = raolal volume, 
ft 5 /lb-mole] 

(also tabulated in 25°F increments over 
temperature range) 


1250-2350 


PVT-closed 
chamber with 
flexible 
diaphragm as 
null detector 


Virial equation of state 
developed using graphical 
evaluation of coefficients 


Ewing et al. 

(1965) 


1 


O 


PV = A(T) + B(T)P + C(T)P 2 + D(T)P : 
A(T) = 0.09268T - 12.4307 

B(T) = 5.80719 - 1 - 8 - 2 - -- + - 68 ^ 2 

t t 3 


C(T) = 2 - 4 ^ 2 ^° + °-30167 _ 1.35708 
t 2 t 4 t 6 


D ( T ) = 0.02146 _ 0.1-6282 + 0.27364 
t 2 t 4 t 6 

[T = °K; t = 10” 3 T;P = atm; V = liters/g-mole] 
(also tabulated) 


1250-2150 


PVT-closed 
cylinder with 
movable probe 
to indicate 
diaphragm 
displacement 


Empirical equation assuming 
PV approaches finite value 
as P ■» 0 


Achener et al. 


IT 


See Table of References. 
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Compressibility Factor 
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Fig. 6. Cesium: Compressibility, Saturated Vapor. 
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ORNL DWG. 67-6492 

















































































































































































































































































































































































































Table 5. 


Enthalpy (Saturated Liquid) 


Equation 

! Eq- 

No. 

1 

Temperature 
Range (°F) 

Standard 

Deviation 

Method 

- [ 

Comments 

Investigator 

(Date) 

* 

Refs. 

Hj, - H 29 g = -1506 + 7.25 T 

[H^p = cal/g-atom; T = °K] 

1 

( 1 ) 

| 

: 

600-1770 

± 1 . 5 1 

Calorimeter 

Not stated whether corrected 
for vapor condensation; 
total data range =170-1770°F, 
below 600 °F an anomoly was 
observed. 

Tepper ( 1963 ) 

12 

H^j, - E 00 = 2.711 + O. 085 4: T - 4.803 

( 2 ) 

152-1656 

0.77 Btu/lb 

m 

Calorimeter 

Corrected for vaporization 

Achener (1964) 

2 

X 10 " 5 T 2 - 1.995 x 10" a T 3 

1 



1 




[Hj = Etu/lb m ; T = *F] 

. ! 

. 







— 

Hj - H d = 3.55 + O.C55 1 * T 

,3, 

: Si 

86-212 

~7 1 

(estimate) 

Calorimeter 

Corrected for vaporization; 
data in range 212-572°F are 
anomolous. 

1 

Lemmon et al. 
(1964) 

4 

Hj - H Q = 0.5 + 0.06 T 


572-2100 






[Hj. = 3tu/lb m ; T = °R] 







1 

Hj - H 0 = - 2.6965 + 0.05683 T 

[1% = Etu/lb m ; T = °R] 

w 

' 

! 

. 

>570 



Based on a value of C which 
is consistent with ^ 3 rd-iaw 
calculations of NRL P-V-T 
dat a. 

Ewing et al. 

( 1965 ) 

1 

Hj - H 2qQ = -1704. 5 - 7. 455 T 

[1^ = cal/mole; T = °K] 

( 5 ) 

l 

1 

1 

600-1770 



Corrected equation based on 
measurements of Ref. 12. 

Heimel ( 1965 ) 

15 


* 

See Table of References. 
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Fig. 8. Cesium: Enthalpy, Saturated Vapor. 























































































































































































































































































































































































































Latent Heat of Vaporization (Btu/lb ) 
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Fig, 9. Cesium: Latent Heat of Vaporization. 


















































































































































































































































































































































































































































































































































































































































































Table 6. Latent Heat of Vaporization 


Equation 

Eq. 

No. 

Temperature 
Range (°F) 

Standard 

Deviation 

Method 

Comments 

Investigator 

(Date) 

* 

Refs. 

h = 257-4 - 0.0362 T 

[h fg = Btu/lb m ; T = °F] 

( 1 ) 


0.48$ 

Direct measure¬ 
ment of heat 
required to 
vaporize a given 
mass of liquid 

• 

Achener (1964) 

2 

h fg = J P s ( 0-53290 Vv g ^ ) 

[h = Btu/lb ; P = atm; T = 3 R; v = ft 3 /lb ] 
fg m s 9 ' n J 

(Tabulated in Appendix A, Ref. l) 

( 2 ) 



Clausius- 
Clapeyron 
equation cal¬ 
culation 

dP 

~ and v , v. from NRL 
dT g’ l 

measurements 

Ewing et al. 

( 1965 ) 

1 

h fg = IT.- 680 . 3 - 2 , 699 .4 ( 10 * 0 ^ 

[h f = cal/g-mole; T = °K] 

( 3 ) 



Clausius- 
Clapeyron 

Investigator’s I-V-T data 

Tepper ( 1965 ) 

17 

= 2.303 BP ( y) 

fg T M g V 

[h^. = cal/g-mole; P = atm: T = °K: 

K = 24.218 cal/liter*atm; 

V = liters; B = 2nd virial coeff] 

(4) 



Clausius- 
Clapeyron 

Second virial coefficient 
tabulated in reference 

Dillon (1965) 

16 


* 

See Table of References. 
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Fig. 10. Cesium: Liquid Entropy 




















































































































































































































































































































































































































































































Table 7* 


Entropy (Saturated Liquid) 


E qua“ ion 

• • 

a* 0 

W k< 

Temperature 
Range (°F) 

Standard 

Deviation 

Method 

Comments 

Investigator 

(Date) 

Refs? 

S - Se 3 7 = -0.8093 + 0.3275 logio T 

-1.5107 X 10“ 4 T + 2.9925 X 10" 8 T 2 

[T = °R; S = Btu/lb ’°R] 

( 1 ) 

152-1656 


• 

r T m dT AH 

- Jb 3 7 p,i T "T m 

+ r T dT 

“T m 3 P,g T 

Achener (1964) 

2 

S - S Q = 0.05683 en T - 0.19387 

[T = °R; S = Btu/lb *°F] 

(2) 

570-2400 



Based or. c = O.O 5683 which 

Pj ^ 

gave gcod degree of internal 
consistency in P-V-T data; 
error in P-V-T iata would 
affect c . and hence S. 

p,£ 

Ewing et al. 

( 1965 ) 

1 


* 

See Table of References. 
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Fig. 11. Cesium: Entropy, Saturated Vapor 
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Fig. 12. Cesium: Liquid 


c .. = 0.05683 Btu/lb- °F 

p^liq m 


For T > 1030°R 


Specific Heat 


[Eq. (5) ] 


< 


* 


Table S. Specific Heat (Liquid) 


Equation 

Eq. 

No. 

Temperature 
Range (°F) 

Standard 

Deviation 

Method 

Comments 

Investigator 

(Date) 

* 

Refs. 

= 7.25 cal/°K*g-atom 

u: 

>600 


Derivative of eq. 
for saturated- 
liquid enthalpy 

No V ^ correction 

Tepper ( 1963 ) 

12 

c = 0.08543 - 9 .6C5 X 10 ' 6 T 

P 

+ 5.985 X 10 - 8 T 2 

[c = Btu/lb ; T = °F] 

L p m’ J 

( 2 ) 

152-1656 


Derivative of eq. 
for saturated- 
liquid enthalpy 

No V ^ correction 
dT 

Achener '1964) 

2 

c = 33-6 - 0.019 T 

P 

[c = joule/g-atom*°C; T = °C] 

( 3 ) 

mp—212 

• 



Rengade (1913) 

18 

c = 0.0554 cal/g*°C 

P 

c = 0.060 cal/g*°C 

P 

00 

86-212 

572-2100 



dP 

No V ^ correction; was 
as much as 2% in places 

Lemmon et al. 
(1964) 

4 

c = C.05683 Btu/lb • °F 

p 31 

(5) 

>570 



Gave greatest degree of 
internal consistency 
in 3rd-law calcula¬ 
tions (trial and 
error) 

Ewing et al. 
(1965) 

1 

= 7.455 cal/mole 1 °K 

(6) 

>600 

* 



Heimel ( 1965 ) 

15 


See Table of References. 
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Fig. 13. Cesium: Vapor Specific Heat. 
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Table 9. Specific Heat (Vapor) 


Equation 

Eq. 

No. 

Temperature 
Range (°F) 

Standard 

Deviation 

Method 

Comments 

Investigator 

(Date) 

* 

Refs. 

( \ r f 

(c p } s - L V s? / p J 3 

[c = 3tu/lb • °F; Ah = enthalpy change 

1? 

(Etu/lb m ); supscript s = saturation] 
(Tabulated in Appendix B, Ref. 1) 


1250-2550 


Calculation 

Normal evaluation of the 
differential at 50°F 
increments 

Ewing et al. 

(1965) 

1 

Superheated Vapor 

Equation giver, in Ref. 1, p. 21; also 
tabulated in Appendix B. 


1250-2550 


Calculation 

Computed along constant 
temperature lines with 
each saturation state 
as starting point. 

Ewing ez al. 

(1965) 

1 


4 

See Table of References. 
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Fig. 14. Cesium: Liquid Surface Tension. 



































































































































































































































































































































































































































E quation j 

Eq. 
No. . 

Temperature 
Range (°F; 

Standard 

Deviation 

Method 

Comments 

Investigator 

(•Date) 

* 

Refs. 

a = 76.*4 - 0.030 (T - 83 ) 

[ct = dyne/cm; T = °F] 


150-1850 

• k% • '* 

Maximum 

bubble pressure 

• 

Tepper et al. 
(1964) 

\ * 

6 , IT 

♦ ♦ 


\J1 
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TRANSPORT PROPERTIES 

Introduction 


The transport properties characterize a material with respect to the 
molecular processes of momentum, heat, and mass exchange. Viscosity, 
thermal conductivity, thermal diffusivity, thermal diffusion ratio, and 
mass diffusivity fall within this category. The electrical resistivity 
is also included in this section, since measurements of this property are 
often used, through the Wiedemann-Franz-Lorenz relationship, to calculate 
the thermal conductivity. 


Viscosity. 


The viscosity, n, characterizes the flow resistance of simple fluids 
and appears as the proportionality coefficient in Newton's law of viscosity: 


dv 


T = —U 

yx 


x 


dy 


. ,■* 

which states that the, shear force per unit area is proportional to the 
negative of the local velocity gradient normal to the mean flow. The 
negative sign specifies that the viscous (molecular) momentum flux tends 
to go in the direction of decreasing velocity. This also enables a treat¬ 
ment which is consistent with that commonly employed for energy (heat) and 
mass transport. 

Liquid Viscosity . The recommended equation [Fig. 15-and Table 11, 

Eq. (l)], obtained by Achener and Boyer, 3 is a reasonable fit to most of 
the available data. The results of Lemmon et al. 4 have been disregarded 
at this time, since the data obtained by these investigators differed 
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significantly from those reported in other sources. . The equation of 
Tepper et al. 1T yields values about 10 $ higher than those predicted by 
the recommended equation. 

Vapor Viscosity . Values along the saturation line are plotted in 
Fig. 1 6; the available data are summarized in Table 12 . The upper curve 
gives the calculated results of Weatherford; 29 the lower curve is the 
polynomial fit [Eq. ( 2 b)] to the data, of Lee and Bonilla. 23 The data of 
the latter are recommended. Note, .however, that the power* [Eq. ( 2 a)] and 
polynomial [Eq. ( 2 b.)] forms given by Lee and Bonilla give somewhat dif- 

. i 

ferent results. Thus: 


T 

(°R) 

Power Fit Polynomial Fit 

(lb m /hr-ft) 

1710 

0.0553 

0.0533 

1800 

O.O568 

0.0568 

1900 

0.0584 

0.0597 

2000 

0.0600 

0,0618; 

2100 

0.0616 

0.0628 

2200 

0.0631 

0.0645 


At this time, no-preference is indicated between' the two forms of the • 
Lee and Bonilla correlation.' Reasonable agreement exists between the 
Kitrilakic and Meeker 19 results, calculated using experimental thermal 
conductivity values and the kinetic theory relationship between \i and k, 
and the Lee and Bonilla data,. On the other hand, the experimental data 
of Weaver and Stratton 22 show significant deviation from that of Lee and 
Bonilla, possibly due to their use of the Taylor and Langmuir 14 relation¬ 
ship for the vapor pressure. 
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Thermal Conductivity 


The thermal conductivity, k, of a material is the proportionality 
coefficient appearing in Fourier's law of heat conduction: 


dT 



This is the one-dimensional form of the equation relating the heat flux 
to the temperature gradient. For an isotropic medium (in which k has the 

same value in all directions) and T = T(x, y, z), the same equation applies, 

♦ 

with the appropriate temperature gradient, in each of the three Cartesian 
directions. In a more general form, the law of heat conduction can be 
written: 

B 3 ! B 2 ! 3 2 T 
3x 2 By 2 3x 2 




neglecting heat sources, where t = time and a = thermal diffusivity 
(k/p <s ). This latter quantity' is often of importance and can be calcu- 

' ST 

lated from the given values for k, p, and c . 

J?. 

Liquid Conductivity . The data of Lemmon et ai . 4 constitute the only 
thermal conductivity measurements for liquid cesium'and are currently 


recommended [Fig. 17 and Table 13, Eq. (3)] the best available values. 
The other tabulated'equations (Tepper et al., 5 Hochman and Bonilla, 18 and 
Hyman 25 ) were obtained through calculation from electrical resistivity : 
measurements; the agreement with the Lemmon et al. results is fair. 

An expression suggested by Ewing and co-workers 30 - for the calculation 
of thermal conductivity from electrical conductivity was applied to the 
resistivity data. However, this procedure yielded results which differed 
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more from the experimental data than the values calculated by the more 
normal Wiedemann-Franz-Lorenz relation. 

Vapor Conductivity , The experimental results of Kitrilakis and 
Meeker 19 [Fig. 18, lower curve, and Table l4, Eq. (l)] are recommended 
on the basis that vapor viscosity values predicted from their data agree 
reasonably with.the best available experimental data on vapor viscosity. 

The calculated values of Weatherford 26 [Fig. 18, upper curve, and Table 14, 
Eq. (3)3 derive from his vapor viscosity prediction assuming a constant 
Prandtl modulus of 0.73-and a frozen specific heat as estimated by Shapiro 
and Meisl. 31 Note that .Weatherford predicts conductivities which are 
about twice those measured by Kitrilakis and Meeker. In contrast, Gottlieb 
and Zollweg 20 reported conductivities measured by a heat-transfer technique 
which are about an order of magnitude smaller than the Weatherford pre¬ 
diction. Achener is currently also using a heat-transfer method to measure 
cesium-vapor conductivities; data have not yet been reported. Cooke (ORNL) 
has proposed a .measurement■utilizing a steady-state variable-gap conduc¬ 
tivity apparatus (an improved version of the technique used by Kitrilakis 
and Meeker) which has been shown to give excellent data for argon and 
helium. 

Electrical Resistivity 

» 

While the electrical resistivity of a liquid is needed in some appli¬ 
cations (e.g., in the design of control and measurement instrumentation)^ 
major interest in this property derives from the relative ease in per¬ 
forming electrical resistivity measurements as opposed to thermal con¬ 
ductivity measurements . The Wiedemann-Franz-Lorenz (W-F-L) relationship, 



4a 

k/p e /T = L o , 

is then used to calculate k from the resistivity, p^; T is the absolute 
temperarature and L q , the so-called Lorenz number* Sommerfeld deduced 
through quantum mechanics that 

< i *■ 

L = tt 2 K 2 /3e 2 = 2.57 X 10“ 8 Btu«ohm/hr* (°R) 2 , 

« 4 * 

where k = the Boltzmann constant and- e = the elementary charge of an 
electron. For molten metals, the W-F-L relation should give thermal 
conductivities close to actual values. However, comparisons made with 
several liquid metals whose thermal conductivities and electrical resis¬ 
tivities are well-known give k predictions higher than the experimental 
values . . It has not yet been resolved whether the discrepancies result 
from a breakdown in the theory or from inaccuracies in the experimental 
measurements. 

Liquid Resistivity . No recommendation as to the best equation is 
made at this time; further evaluation of data in accord with the above 
discussion remains to be carried out. k 

The equations of Hochman and Bonilla 18 and of Kapelner and Bratton 24 
agree closely at 1700°F (~198 pohm*cm versus ~19^ pohm-mc). The results 
of Tepper et al. 5 fall about k°/o below these (~l88 [lohnrcm) and of Hyman 25 
about 28$ low (~1^2 jaohm*cm). Either of the first two equations [Eqs. (2) 
and ( 3 ) in Table . 15 ] appears adequate; the results of Hochman and Bonilla 


are shown in Fig. 19- 



Mass■Diffusivity 


While no data are included in this collection at the present time 
on the mass diffusivity, D, for the sake of completeness a few brief 
comments are merited. In a fashion analogous to that for viscous flow 
or heat conduction, the diffusivity is the proportionality coefficient in 
Fick's law of diffusion which, in one-dimensional form, can he written: 


dC 


j Ay D AB 


A 


dy 


where is the molar diffusion flux of component A in the y-direction, 

C. is the molar concentration of component A in the mixture, and D (= D ) 

■H. n KA. 


is the diffusivity coefficient for component A through component B. Of the 
three transport processes, diffusion is the most complicated since it deals 
with mixtures. . This property assumes importance in respect to the compati- 
bility of cesium with its container materials, condensation of cesium vapor, 
etc. 


Thermal Diffusion Ratio 

This describes the tendency for species to diffuse under the influence 
of a temperature gradient. Generally, this effect is quite small; but 
systems have been developed which produce temperature gradients sufficiently 
steep to effect separation of mixtures. The thermal diffusion ratio is 
defined as: 



where p is the fluid density; 


c 


is the total molar concentration; 



and 




42 


T 

are the molecular weights of species A and B, respectively; is the 
thermal diffusion coefficient for species A in the mixture; and is 
the binary diffusivity in the mixture A-B. A discussion of this property 
and the interrelationships between the transport properties can be found 
in Section 18.4 of Ref. 32* 
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Fig. 15. Cesium: Liquid Viscosity. 


















































































































































































































































































































































































































































































































































































































Table 11. Viscosity (Liquid) 


Equation 

Eq. 

No. 

Temperature 
Range (°F) 

Standard 

Deviation 

Method 

Comments 

Investigator 

(Date) 

* 

Refs. 

logic U = l8 y 23 - 0.59911 logio T + I.U1U67 

(1) 

mp-1958 

3.516 

Oscillating 

Cylinder 

Used density of Ref. 2 

Achener and 
Boyer (1966) 

3 

Cki = lb m /hr- ft; T = °R] 








logio H - 0.401 

(2) 

85-620 


Oscillating 

Cylinder 

Used density of Ref. 2 

Lemmon et al. 
(1964) 

4 

logio H 


620-2100 






[p. = lb /hr*ft; T = °Rl 
m 




. 

%k 



„ = 10.14 e 576/r 

(3) 

83-412 


Oscillating 

Sphere 


Andrade and 
Dobbs (1951) 

21 

[p = millipoise; T = °C) 








logio u = —j 5 '- - 0.641-91 

(4) 

99-2000 


Oscillating 

Cylinder 


Tepper et al. 

(1965) 

IT 

[|i = lb m /hr-ft; T = °R] 









* 


See Table of References. 
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Fig. 16. Cesium: 


Vapor Viscosity. 






























































































































































































































































































































































































































































































































































































Table 12. Viscosity (Vapor) 


Equation 

Eq. 

No. 

Temperature 
Range (°F; 

Standard 

Deviation 

Method 

Comments 

Investigator 

(Date) 

Refs.* 

|i = 0.99 X 10" 9 T 1,73 

[|i = go/cm- sec; T = °K 

(1) 

98O-31IK) 


Transpiration 


Weaver and 
Stratton 

(1965) 

22 

1 i = 6.2887 T 0,6241 

or 

(2a) 

Il60-l8l0 

± 3-996 

Coiled Capillary, 
Steady Flow 


Lee and 

Bonilla 

(1966) 

23 

“ ■ ^ 48 - 14 + 1511 - 52 ( iroo ) 



±3.2* 





^ 28 - 65 ( 1000 ) 

(2b) 







[p = micropoise; T = °K] 








p = 7-05 W 2 

[p = micropoise; T = 3 K] 

(3) 



Derived 

Based on kinetic theory 
relationship between p 
and k; calculated from 
k measurements 

Kitrilakis 
and Meeker 

(1963) 

19 

Curve plotted in Ref. 26. 

* 

W 

5*+0—22*40 


Calculated 

Computed for variable 
molecular weight 
saturated vapor using 
new equation for gas 
viscosity 

Weatherford 

(1961) 

29 


See Table of References. 
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Thermal Conductivity (Btu/h 
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Fig. 17. Cesium: Liquid Thermal Conductivity. 

































































































































































































































































































































































































































































Table 13. Thermal Conductivity (Liquid) 


Equation 

Eq. 

No. 

Temperature 
Range (°F) 

Standard 

Deviation 

Method 

Comments 

Investigator 

(Date) 

Refs.* 

k = 0.273508 + 7.3lK)894 X 10" 6 C 

- 7.422601 x 10" 8 T 2 + 7.491586 

X 1C" 12 T 3 

[k = w/in.*°F; T = °F] 

:d 

168-1945 


Derived 

Calculated from electrical 
resistivity measurements 
using Weidemarn-Franz- 
Lorenz equation 

Tepper et al. 
(1964) 

5 

k = 0.124664 - 2.51402 x 10’ 5 T 
- I.63806 x 10" 9 T 2 
[k = w/cm*°F; T = 3 F] 

(2) 

600-3000 

2.4# 

Derived 

As above; W-F-L constant = 

2.3 X 10 s (volts/°K) 2 

Hochman and 
Bonilla (1965) 

18 

k = 16. 220 - 0.003995 T 

[k = Btu/hr* ft* °F; T = °R] 

(3) 

347-2160 

• 

Steady-state 
longitudinal 
heat flow 
apparatus 

Least-square fit of data 

Lemmon et al. 
(1964) 

4 

k = O.C48 to 0.050 cal/sec*°C*cm 

w 

mp—786 


Derived 

Calculated from electrical 
resistivity data through 
W-F-L relationship 

Hyman (1961) 

25 


* 

See Table of References. 
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Fig. 18. Cesium: Vapor Thermal Conductivity. 
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Table l4. Thermal Conductivity (Vapor) 


Equation 

Ec. 

No. 

Temperature 
Range (°F) 

Standard 

Deviation 

Method 

Comments 

Investigator 

(Date) 

* 

Refs. 

k = 1.65 X 1CT 6 T 1 / 53 

g 

[k = w/cm*°K; T = °X] 

a) 

931-3184 


Steady-state, 
parallel plates 


Kitrilakis and 
Meeker (1963) 

19 

T(°K) k(w/cm*°K) 

945 1.l4 x 10" 5 

1000 1. i-9 

1060 l. 75 

1123 2.C 

1188 2.2 

(2) 

1241-1679 

. * 

Heat-loss 

Method of calculation not 
clear 

Gottlieb and 
Zollweg (1983) 

20 

Curve shown in reference 

(3) 

54o-224o 


Calculated 

Obtained using frozen 
specific heat and 
viscosity (Ref. 29) 
assuming a constant 

Pranitl modulus of 

0.73 

Weatherford 
et al. 

26 


* 

See Table of References. 
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Electrical Resistivity (p.ohm* cm. 
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Fig. 19. Cesium: Electrical Resistivity. 

































































































































































































































































































































































































































Table 15. Electrical. Resistivity (Liquid) 


Equation 

Eq. 

No. 

Temperature 
Range (°F) 

Standard 

Deviation 

Method 

Comments 

Investigator 

(Date) 

Refs.* 

p = 10.90857 + 3.39024 x 10' 2 T - 1.670142 

v» 

X 10“ s T 2 + 1.096439 X 10" 8 T 3 

[T = °F; p = pohm*in. ] 

6 

(1) 

168-1946 

1% 


Measured resistance and 
calculated resistivity 

Tepper et al. 
(1964) 

5 

in = 3.4630 + 1.4872 X 10" 3 T 

- 3.62046 x 10’ 7 T 2 + 6.99129 

X 10- 11 T 3 

[T = °F; p^ = pohm* cm] 

(2) 

600-3000 

2. % 

• 

Measured conductance and 
calculated resistivity 

Hochman and 
Bonilla (1964) 

18 

P. = 34.383 + 1.124 x 10- 1 T - 2.869 X 10‘ E T 2 

• V 

+ 1.007 x 10~ 7 T 3 

[T = °C; p = pohm*cm] 

0 

(3) 

83-16IC 

±0.24 % 


Measured potential drop 
with known current 

Kapelner and 
Bratton (1962) 

24 

P e = 38.5 [1 + 0.00308 (T - 50)] 

[T = °C; p e = piohm- cm] 

w 

mp-786 

■** '•«» 

±1% 


% 

Hyman (1961) 

25 


* 

See Table of References. 
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A, B, C, D, E 


NOMENCLATURE 

Constants appearing in equations of state (Table 4); 
also designates species in a mixture 

Total molar concentration; C., molar concentration of 
species A in a mixture 

Heat capacity; c p ^ g , of vapor; c p ^ of liquid; (c ) g , 
at saturation conditions 


Coefficient of mass diffusivity; D._. (= D_. ), binary 

AB BA 

diffusivity in a mixture containing species A and B 
Thermal diffusion coefficient for species A in a mixture 
Enthalpy; H_, at temperature T; H , at 0°R; H , at 

•rjTi 4 O 29B 

298 °K; H , of vapor at temperature T; H , of vapor 
S E 

at 0°R 


Ay 


P, P 


Latent heat of vaporation 

Molar diffusion flux of species A in y-direction 

Thermal conductivity 

Thermal diffusion ratio (p. hi) 

Lorenz number (= k p /T) 

o 

Molecular weight; M^ of monomer; M ., M^, of species A 
and B, respectively, in a mixture 

Pressure; P 0 at critical condition 

Heat flux in y-direction 


Gas constant 


T o 

Entropy; S , of vapor at temperature T; 8 , of vapor 

at 0°R 8 S 


Temperature; T c , at critical condition; £T , temperature 
difference 


v, V 


Volume; V c , critical condition; v , of vapor; v^, 
of liquid; V, molal volume ® 

Velocity in x-direction 


i 
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x > y, z 

z 

a 

€ 

K 

M* 

P 



Cartesian coordinate directions 
Compressibility factor. (= p V/RT) 

Thermal diffusivity (= k/p c ) 

Jr 

Elementary charge of an electron 

Boltzmann constant 

Viscosity 

Density 

Electrical resistivity 
•Surface tension 
Time 

Shear stress^ or viscous flux of x-momentum in y-direction 
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